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Abstract. Cellular automata can simulate many complex physical phenomena using the 
power of simple rules. The presented methodological platform expresses the concept of pro- 
grammable matter in which Newton’s laws of motion are one of examples. Energy has been 
introduced as the equivalent of the ” Game of Life” mass, which can be treated as first level of 
approximation. The temperature presence and propagation was calculated for various lattice 
topology and boundary conditions by using the Shannon entropy measure. The conducted 
study provides strong evidence that despite not fulfillment the principle of mass and energy 
conservation, the entropy, mass distribution and temperatures approaches thermodynamic 
equilibrium. In addition, the described cellular automata system transits from positive to a 
negative temperatures that stabilizes and can be treated as a signature of system dynami- 
cal equilibrium. Furthermore, the system dynamics was presented in case of few species of 
cellular automata competing for maximum presence on given lattice with different boundary 
conditions. 


1 Introduction to Classical Conway’s Game of Life (CCGoL) 


A cellular automaton is a system consisting of cells arranged most often on a one, two or three 
dimensional regular lattice, which at a given moment are in one of M-th possible states expressing 
M-valued logic. The dynamics of the model depends on the definition of individual cell states and 
the rules of transitions between them [13]. One of the simplest examples is a one-dimensional cellular 
automaton. Suppose that the cells placed on the lattice can be in one of two states, which are marked 
with white (default assigned to dead state or logical zero) or black color (default assigned to alive 
state or logical one). We define a rule that if a given cell is black, then the cell to the right of it will 
change its state. This situation is depicted in a Figure|1} We can observe how the system changes 
in the subsequent steps of the simulation. The parameter determining the change of the cell state 
is the state of the left neighbour of a given cell. There are many other possibilities to choose such a 
parameter, e.g. the condition of the state of both neighboring cells or having the nearest neighbors 
with opposite states. In order to determine system dynamics, we must have a defined initial cellula 
automata states (information about initial system dynamical state) and a specific set of deterministic 
or probabilistic rules. 
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Fig. 1: Evolution of a one-dimensional cellular automaton in successive cycles with left side partial 
logical negation rule (if the state of nearest left cell is alive then given automata cell change its state 
to opposite). 


The Conway’s Game of Life [3] is an example of a cellular automaton with a deterministic rules. It 
was proposed in 1970 by John Conway and cellular automaton system consists of cells located on 
a two-dimensional lattice, which can be in one of two states: alive or dead. The rules specify the 
required number of neighbors and cell states that are taken into account to determine their state in 
the next cycle. Given the neighborhood through the 8 closest cells, we can write the 3 main rules of 
the Classical Conway’s Game of Life (CCGoL): 


1. If a dead cell has exactly 3 neighbors, it comes alive in the next cycle. 
2. If a living cell has 2 or 3 neighbors, it survives in the next cycle. 
3. If a cell has a different number of neighbors than stated above, it will be dead in the next cycle. 


The rules defined in this way allow for the generation of various types of structure topologies with 
automaton state set to 1, as shown in Figure [2] The most common are ”unstable structures”, which 
change in successive cycles, but do not return to their initial state. A single cell cannot survive on the 
lattice, because it has less than 2 neighbors. Dead cell surrounded by alife cells cannot come to life, 
because it has number of neighbors is different from 3. If the simulation is continuing sufficiently 
long time, on the lattice usually remain structures that are unchanging in time - "still lifes” (an 
example is the ” block” shown in the Figure P) or changing in time in periodic way, so they return to 
their original shape after k cycles - ” oscillators” (an example is the ” blinker” shown in the Figure[2}. 
There are also structures that move in a certain direction - ” gliders” (Figure [3) that are analogical 
to 1st Newton dynamics preserving momentum (speed and direction of propagation in this case), 
leave a trace of ” blinkers” - ”star ships” and objects, which periodically eject ” gliders” - ” guns” [8]. 
Figure [4|shows one among many existing in the Conway’s Game of Life oscillators during successive 
iterations of the system simulation. 

The ”toad” oscillator has a period of 2, which means that it has continuous switching between 2 
different fixed configurations. Each 2-dimensional discrete lattice field has a specified number that 
indicates the number of neighbors of the given cell. If the number is green, the cell will be alive in 
the next cycle. If the number is blue, that cell will be dead in the next cycle. 
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Fig. 2: Evolution of various topologies of cellular automata structures in time with deterministic 
rules of CCGoL. One can identify two dynamically unstable structures and two structures that have 
dynamical stability with time. 


(a) (b) (c) (d) (e) 
Fig. 3: Evolution of the ” glider” configuration of cellular automata having propagation property with 
time in deterministic CCGoL. 


Fig. 4: Evolution of the toad” configuration of cellular automata being an oscillator in successive 
cycles in CCGoL. 


2 Introduction to Stochastic Classical Conway’s Game of Life (SCCGoL) 


The Stochastic Classical Conway’s Game of Life (SCCGoL) was created by an addition of a cell spon- 
taneous change probability to the rules that were initially deterministic, so stochastic determinism 
was achieved. With a given prefixed spontaneous probability value, the state of the cell can change 
regardless of the number of neighbors. Cell states have values between 0 and 1 and are called mass. 
Due to the fact that SCCGoL have different rules from CCGoL, cells have almost never exactly two 
or three neighbors. A condition for given cell to come to life from a dead cell state (creationism of 
alife cell) is by having a number of neighbors in a certain range of values. Similar rules applies to 
a living cell justifying its alife or dead state in next time iteration. By setting standard intervals of 
allowed/forbidden number of neighbours values, in which cell is alife/dead and by adding the addi- 
tional spontaneous rule probability for cell to change in the next iteration (probability of changing 
the state of a cell regardless of the number of neighbors), we are able to create a simulator, where the 
cells never die or it is very difficult from a probabilistic point of view for a given cell to stay alive. 
If, maintaining standard neighbours condition for being dead/alife in previously defined intervals 
and having a selected initial cell alife automata configuration, we can systematically increase the 
spontaneous probability level from 0% to 100%, and we result in the graph as depicted in Figure 
Simulations were conducted for a lattice size 10 by 10 with the initial condition of a cellular 
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Fig. 5: (a) Schematic view of initial conditions in SCCGoL. (b-c) Dependence of automata population 
average cycle lifetime (over 1000 trials) on probability of spontaneous change of cell state from life 
to dead and reversely with preservation of standard rules in Conway’s Game of Life. 


automaton 2 by 2 (Figure |5h) with limited maximum number of cycles set to 1000 in conducted 
simulations. For a probability of changing the state of a cell regardless of the number of neighbors 
equal to zero, there is full determinism, therefore it is a situation of CCGoL with changed numbers of 
neighbors. As the probability increases (in range from 0 to 2%), the life expectancy of the population 
decreases due to too few neighbors. Starting simulation from a probability level close to two percent 


(probability of changing the state of a cell regardless of the number of neighbors), the average life 
expectancy of the population begins to rise, which is caused by the more frequent appearances of 
living cells. Conducting the simulation with a probability level above eighteen percent we observe 
that the population practically never dies and keeps average cycle life time being at least 1000 or 
more time iterations. 


3 Generalization of Stochastic Classical Conway’s Game of Life 
(SCCGoL) to the case of N competing cellular automata species 


The created simulation platform enables to create N different cellular automata species that com- 
petes between themselves for the resources enabling them to replicate. The given automata species 
has better replication properties within its own community and worse replication properties in case 
of neighbors being of different species. Such situation is normally encountered in human population 
when people of one homogenous identity prefer to collaborate more than in the case of people hav- 
ing much different identity (being from different tribes). Due to that fact soft antagonistic relation 
between different species is introduced indirectly by means of higher level of tolerance (or effective- 
ness of replication) towards neighbors of own species than neighbors of other species (other different 
species are treated in the same way). In real way this is an analogy of members collaboration of a 
given nation or culture within given culture or community versus different culture or community. 
Let us consider N = 4 (number of different automata species), so we have the following determined 
formula for number of effective existing neighbors Neef fective,k for given k-th tribes as: 


Neef fective,k = aik Ne + a2 k Nez +... + ak,kNek +... + ak, NNen. (1) 


Previously defined replication rules promoting its own species can be formally expressed by following 
condition (max(as k) = ap k and akk > asp if s # k). In the conducted simulations all tribes have 
assigned the same value (41,1 = d2,2 = 43,3 = 44,4 3) and another same value for (a1,2 = 41,3 = 
a1,4 a2,1 a2,3 a2,4 43,1 a3,2 a3,4 a4,1 a4,2 a4,3 3) what simply means 
that automata tribes promotes its own tribe and only partly promotes other tribes with no special 
distinction on which different tribe it is pointing to. Cellular automata species distribution for k-th 
tribe across two-dimensional lattice is generated with use of a two-dimensional Gaussian distribution 
given as follows: 
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where f(x,y), is a mass function depending on the cell coordinates, (£o,k; Yo,k) the center of the 
given cellular automaton species. Given parameters Ag, and Or k, Cy,k can control the mass of 
species and spread across the ”Globe” as depicted in Figure a In the case of several species on 
the lattice, SCCGoL simulation results promote the formation of new cells among the species that 
has the most mass in the neighborhood as in accordance in the Figure [7] During life cycle always 
newly formed cell has a mass that is more dependent on cells of the same species than other species 
of its neighbors. We observed that automata species fight each other in order to spread across the 
lattice and exclude other species what is an expression of some form of biological Darwinism being 
a consequence of formula [I] 


(a) 
Fig. 6: (a) Case of mass distribution of cellular automata of one species using a two-dimensional 
Gaussian function being isotropic. (b) Example of initial distribution of three cellular automata tribes 
with the case of boundary existence between each separate tribe with a rule that one geometrical 
place on the lattice is occupied by cellular automata of given species with dominant mass. 


(a) (b) 
Fig. 7: Evolution of map distribution of 4 cellular automata populations of different species with 100 
by 100 lattice size. One can spot effective tendency of each species to occupy maximum possibly 
territory at the cost of other species in 4-species SCCGoL (SCCGoL4S) what can be understand as 
occurrence of weak antagonistic relation 


4 Methodology of description of SCCGoL dynamics by tools of classical 
statistical physics 


Adding probability to the initially deterministic Game of Life and changing the interaction rules 
between neighboring automata brought the necessity for a new quantity called mass that has con- 
tinuous real values (other than 0 and 1 present in deterministic Game of Life). At this stage vividness 
of the whole cellular automata population can be understood and approximated by the population 
total energy (where simply mass is equal to energy). The level of the automata distribution order is 
expressed by the entropy. Entropy was introduced by Rudolf Clausius in 1865 as a thermodynamic 
state function [II]. If the entropy at the initial state and the entropy of the final state are devoted 
by S; and Sy respectively, then we have: 


fdE dS 1 
s-3= | T dE T (3) 


Very last formula gives definition of a temperature under assumption that energy and entropy is 
defined that is the case. We use instead of thermodynamic entropy the Shannon entropy measure 
given as: 


Sshannon = X v(ai, yj) log p(x, Yj) =E |- log p(z, y)| (4) 


ij 


In order to calculate the probability in the equation [4] the SCCGoL simulation is repeated several 
hundred times, so one determines an average population cell mass in successive cycles. The calcula- 
tion of the temperature, which is a measure of thermal state was achieved with an equation [3] with 


a changed form: 
dE(«,y,t) 
eres dE t 
T(a,y,t) = : ae 
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The temperature defined by last formula can be calculated by two different methods. The first 
method relies on the calculation of the energy and entropy for the entire system, and then numerical 
calculation of the derivative of energy in function of entropy. The second calculation approach is 
to differentiate mass and entropy with respect to simulation time for each cell and as a result of 
corresponding ratio we obtain a temperature map. The example of evolution of mass and entropy 
for SCCGoL is depicted in Figure |$] and shows maximization and saturation of entropy, what is 
confirmation of Second Thermodynamic Law that is valid in physical systems, while we are dealing 
with cellular automata system. On another hand the mass of a system saturates and increases to 
certain critical value as given in left part of Figure [5] 
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Fig. 8: Evolution of mass and entropy in successive cycles in SCCGoL manifesting maximization with 
saturation and approaching stationary thermodynamical equilibrium with initial condition depicted 
in Figure pl 


5 Numerical analysis of SCCGoL dynamics with methodology of 
classical statistical physics 


Numerical simulations were carried out for 4 topologies cellular automata (case of Figures (Dh, mp, 
[13h, [15p). We preinpose such a rule that given alive cell has 20% probability of changing its state 
to dead state and that initially dead cell has 20% probability of changing its state to alive with 
a mass choose randomly from an interval value in (0,0.5). Still given cell state is depended on its 
neighbors, since it has 80% probability of changing its state due to state of neighbors. In order 
to obtain cellular automata probability map dynamics, simulations average of cellular automata 
positions was conducted. Figure [10] shows results that were obtained for a lattice size 100 by 100 
with a one cell alive as initial lattice state (case of Figure 9h). With subsequent cycles, the cells 
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(a) 
Fig. 9: Evolution of diffusion process in cellular automata system for a limited lattice size 100 by 100 
in SCCGoL (averaged over 1000 trials). It shows final saturation of mass and entropy (as depicted 
in Figure|g) what implies approaching thermodynamical equilibrium with characteristic fluctuations 
of mass and entropy around effective stationary values. 


occupy more and more space on the lattice, which can be seen as increase the mass of the entire 
system (possible mass creationism is inherent feature of Conway’s Game of Life). The sum of the 
masses of all cells in successive cycles is depicted in Figure |8} with a comparison of the entropy 
change of the entire system. As we observe in Figure [1Op, high entropy occurs at the edges of the 
population, which is caused by the entropy wave that meets area with almost zero cell occurrence. 
Before equilibrium is established, the entropy of the system slightly decreases due to the loss of this 
extra entropy at the edges as can be seen in the right part of Figure [8] Having established those two 
quantities, we conduct their differentiating with respect to time, and by formula|3] one can establish 
the whole effective temperature of system by dividing change of mass at given time by change of 
the entropy. Figure shows the greater susceptibility of the entropy change to the constraints 
associated with the limited size of the simulation lattice that is ended with impenetrable walls. The 
dependence of mass derivative with respect to time simulation (case of Figure w) does not have 
such large oscillations as in the case of entropy derivative with simulation simulation (case of Figure 
). The temperature calculated by such procedure gives values just above zero (slightly positive) 
up to the 50’th cycle. We can see a large down peak caused by a slowdown in entropy increase. 
From about the 75’th cycle, the temperature of the system goes from positive to negative values. 
A Figure [10}) describes anomalous thermalization process in SCCGoL with a case of approaching 
temperature and entropy equilibrium. Surprisingly, the thermodynamic equilibrium is achieved for 
a case for negative temperatures. A second possible approach in determination the temperature is 
by use the derivatives of the mass and entropy with respect to time of the individual cells and by 
obtaining a temperature map of all cells. The temperature depicted in Figure [10} is mostly negative 
and steady, what corresponds to a situation where mass and entropy have come to equilibrium. We 
can distinguish two regions in the simulation, the first one with no temperature gradient and zero 
negative temperature and the second region with non-zero temperature gradient that also include 
positive temperatures. The place, where time derivatives of mass and entropy have noticeable values 
is at the edges of automata population, where we observe slightly positive temperature, which 
corresponds to the situation in Figure [10} before the 75’th cycle. 
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Fig. 10: Dynamics of thermodynamical parameters (mass, entropy and temperature) with simula- 
tion time in SCCGoL (lattice size 100 by 100) also given in Figure [8] Achieved thermodynamical 
equilibrium is accompanied with final distribution of negative temperature (starting from positive 
temperature distribution as in accordance with formula |3) as experimentally observed necessary 
criteria for final thermodynamical stability. One can spot various similarities of statistical behaviour 
of SCCGoL with physical systems described by classical thermodynamics (maximization and satu- 
ration of entropy, decay of temperature gradients and final thermalization, uniform distribution of 
mass and energy). 


6 Numerical study of one species cellular automata with various 
boundary conditions 


Next family of simulations were conducted with use of barriers impenetrable by cellular automata 
cells via which cellular automata cannot interact with each other. We consider a single cellular 
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Fig. 11: Diffusion process in SCCGoLp20L100b100 (lattice size 100 by 100, 20% probability of 
spontaneous change of cell state from life to dead and reversely with preservation of standard rules 
in Conway’s Game of Life as also given in Figure [5) case of system of two weekly interconnected 
chambers by means of two small holes in barrier. Two stages of diffusion can be spotted in mass and 
entropy dynamics that has two consecutive processes: full diffusion in the left chamber is leading to 
full diffusion in the right chamber. Further details of thermodynamical parameters space dependence 
evolution with time are given by Figure [12] 


automata (single automata seed) placed in empty chamber with impenetrable walls with two small 
holes that link it to another empty chamber as depicted in the Figure Mh. We observe a diffusion of 
cellular automata with simulation time that consists of two main processes: creation and diffusion 
of cellular automata in the first chamber and diffusion of cellular automata from the first chamber 
into second chamber accompanied with creation new automata in the second chamber. Those two 
consequent processes are accompanied by effective slowdown in diffusion that is seen in left part 
of Figure m}. At the same time we observe slowdown in entropy increase as in an accordance 
to the right part of Figure m}. Once mass saturation was obtained in the simulation we observe 
maximization and small drop in entropy that later stabilizes and saturates, as in the right part of 
Figure m. Entropy unlike mass is characterized by a large variation of values in the middle of the 
population and at the edges of the cellular automata population. In the situation with barrier (case of 
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Fig. 12: Dynamics of thermodynamical parameters with simulation time in SCCGoLp20L100b100 
with two weekly interconnected chambers by two small holes in barrier that were initially depicted 


in Figure 


Figure [12p) this is particularly evident after the cells pass through the gaps. This particular process 
is due to the logarithm function dependence of Von Neumann entropy, which tends to negative 
infinity for arguments going to zero from the right. Observed criteria of an equilibrium is fact that 
mass and entropy have steady values and that temperature have negative value in thermodynamical 
equilibrium. Almost always before the equilibrium moment is achieved, time derivatives of mass and 
entropy have positive values, and still the temperature is positive. From the moment the equilibrium 
is achieved, we are dealing with small fluctuations of entropy and temperature. We observe the 
correlation stating that slight increase in mass results in slight decrease in entropy and vice versa. In 
an analogical way to the system with one barrier, simulations were conducted on the system with two 
barriers and the initial structure as depicted in Figure [13h. We consider a single cellular automata 
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(c) 
Fig. 13: Dynamics of diffusion process for a system SCCGoLp20L100b100 with two barriers with 
two small holes in each barrier (generalization of situation from Figure [12) creating three weakly 
interconnected chambers perturbed by mutual interactions mediated by holes. Monotonicity in in- 
crease of entropy is twice shortly interrupted by small decline, what is associated with automata cells 
” colliding” with barriers and experiencing short lasting slowing down in its propagation. Details on 
space depended thermodynamical parameter evolution with time are given by Figure [14] 


placed in empty chamber with impenetrable walls with two small holes that link it to the second 
empty chamber, which is connected to a third empty chamber by impenetrable walls with two small 
holes as depicted in the Figure E3. We observe a diffusion of cellular automata with simulation time 
that consists of three main processes: creation and diffusion of cellular automata in the first chamber, 
diffusion of cellular automata from the first chamber into second chamber accompanied with creation 
new automata in the second chamber and diffusion of cellular automata from the second chamber 
into third chamber accompanied with creation new automata in the third chamber. Twice cells try 


to reach impenetrable barriers, we observe small drops in entropy. After the second time entropy 
stabilizes and saturates, as in the right part of Figure [T3p. In a system with two barriers, it lasts 
longer for mass and entropy to reach equilibrium than in the case of a system with only one barrier. 
As depicted in Figure (14k, due to the cells approaching the barriers and losing the extra entropy at 
the edges of the population, we observe significant fluctuations in the time derivative of the entropy. 
This results in large peaks seen in the Figure E4. 
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Fig. 14: Space depended dynamics of thermodynamical parameters with simulation time in SCC- 
GoLp20L100b100 with three weekly interconnected chambers by four small holes also depicted in 


Figure 
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7 Numerical study of two species cellular automata in perturbative 
interaction by narrow constriction 


Further simulations for the case of two species cellular automata were carried out for a system divided 
into two reservoirs separated by two impenetrable barriers with one small hole (case of Figure [15p) 
that implies perturbative interaction between tribes. In the left upper corner of the first part of 
the system there have been located cells of the first cellular automata tribe. At a closer distance 
to the hole in impenetrable wall, but in the second right reservoir there have been located cells of 
the second cellular automata tribe. As depicted in Figure E5}, we observe similar final masses and 
their dynamics in case of both tribes, but in accordance to Figure (15H, different entropy dynamics. 
Very last is due to the distance of the cellular automata tribes from the small hole in impenetrable 
wall: a tribe located further from that hole needs more time to propagate and occupy its natural 
neighborhood and first left chamber of the system, and thus this tribe has a lower probability of 
taking over the territory of the other tribe. However the noticeable fact is that mass and entropy 
of both tribes achieves equilibrium and finally tribes end up in bit same geometrical and dynamical 
situation. As depicted in Figure [L5F, the right tribe closer to the small hole occupies its nearest 
neighborhood territory more quickly, resulting in an attempt to occupy a rival tribe territory. As 
depicted in Figure [16] we observe large oscillations in the time derivatives of mass and entropy of 
both cellular automata tribes. In contrast to previously conducted simulations, the temperature of 
the system after reaching equilibrium is not only characterized by negative values. Tribes existential 
competition is the reason of occurrence of both positive and negative temperatures. 
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Fig. 15: Diffusion process in case of system with two cellular automata tribes weekly interacting 
with each other via a small hole in double barrier as depicted in (a), where initial configuration 
is presented. After long thermodynamic equilibrium is achieved as given by (b) so two cellular 
automata tribes coexist in two different geometrical domains effectively geographically separated. 
In case of both tribes mass and entropy saturates having tendency to oscillate in thermodynamical 
equilibrium. 


Derivative of mass with respect to time of the first tribe 


Derivative of mass with respect to time of the second tribe 


=z 


400 


200 


—200 


—400 


14 


12 


Fig. 16: Dynamics of thermodynamical variables for the case of two competing cellular automata 
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(c) Temperature with time for first and second automata tribe 


tribes (first tribe is placed on the left and second tribe is placed on the right). 


8 Conclusions and future perspectives 


Cellular automata can simulate many complex physical phenomena using the power of simple rules 
as it was shown in the case of cellular automata diffusion dynamics confirmed by various simulations 
for one and many automata species. Certain type of automata Darwinism was spotted by studying 4 
automata species dynamics as given by Fig[7] The SCCGoL dynamics study provides strong evidence 
that despite the fact that the principle of conservation of mass is not fulfilled, since we have creation- 
ism and annihilation of automata, the entropy and temperature comes to equilibrium. In conducted 
various simulations of Stochastic Conway’s Game of Life dynamics we report transition from positive 
to negative values of temperatures and we are aware that there is maximum level of mass and energy 
density allowed for cellular automata, since otherwise they would die due to overpopulation. The fact 
that the temperature can be negative is known in condensed matter physics, but with assumption 
that the energy is top-bounded. In most ”normal” situations this is impossible, but in a rare cases 
in solid state physics approximately it can be achieved by inverting the population state. Obviously 
there is a such limitation on top-bounded energy value (mass density value) in Stochastic Conway’s 
Game of Life. Therefore, it is still consistent with thermodynamics methodology, as it was pointed 
by professor Adam Bednorz (Faculty of Physics, University of Warsaw). 


Following conclusions were derived basing on conducted simulations and described methodology: 


1. Identification of thermodynamically defined temperature as proper measure of system evolution 


with ’-’ sign (case of Figures[10} (14h. 


2. Identification of mass as effective energy of system (in first approximation) (case of Figures 


MIENE) 


3. Identification of Shannon Entropy as effective system entropy (in first approximation) (case of 
Figures t5}. 

4. Generalization of Stochastic Conway Game of Life of N tribes (approximated analogy to N-body 
Quantum Physics can be conceptionally drawn) as depicted in Figures 

5. Confirmation validity of second law of thermodynamics in SCCGoL (entropy maximises and 
saturates, case of Figures (13). 

6. Identification of short lasting Shannon entropy peak that later minimizes and saturates in SGoL 
(case of Figures|} ma}. Monotonicity in increase of entropy is twice shortly interrupted by small 
decline, what is associated with automata cells colliding” with barriers and experiencing short 
lasting slowing down in its propagation (case of Figure[13). 


Conducted analysis of Stochastic Game of Life allows to treat such system as mathematical object 
well described by methodology of classical statistical physics. Obtained numerical results by various 
simulations suggest that we shall introduce another definition of temperature in Stochastic Conway’s 
Game of Life system by adding ’minus’ sign to temperature known in statistical physics, so we obtain 


the following formula: 
dE 
TemperatureConway—Pomorski—Kotula = TAS (6) 


TemMmperatuTEsStatisticalPhysics =+ 


Having such a definition of Conway-Pomorski-Kotula temperature we can use tools of statistical 
physics in Stochastic Game of Life preserving most classical physics thermodynamical intuition about 
various situations we can come across. There are well-known inherent analogous between classical 
statistical physics [2 [7] [2] [5] [4] and quantum mechanics [I]. Therefore further research perspectives 
in study of Stochastic Classical Conway’s Game of Life assume the usage of quantum mechanics 


being able to simulate classical statistical physics (as expressed by epidemic model or stochastic 
finite-state machine) as explicitly represented by tight-binding model [9] [10] or Schroedinger model 
directly proposing structures implemented in semiconductor single-electron devices. Noticeable var- 
ious obtained map of probability of Stochastic Conway’s Game of Life can be parameterized by 
non-linear Schrodinger equation and especially by Ginzburg-Landau model [6], what can be the 
base for quantization of Stochastic Classical Conway’s Game of Life. 
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